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A Mn-porphyrin was contacted on Au(lll) in a low-temperature scanning tunneling microscope 
(STM). Differential conductance spectra show a zero-bias resonance that is due to an underscreened 
Kondo effect according to many-body calculations. When the Mn center is contacted by the STM 
tip, the spectrum appears to invert along the voltage axis. A drastic change in the electrostatic 
potential of the molecule involving a small geometric relaxation is found to cause this observation. 

PACS numbers: 73.63.-b, 73.63.Rt, 68.37.Ef 


The spin states of molecules at surfaces are intrigu¬ 
ing because of the interaction of localized molecular or¬ 
bitals with delocalized substrate states. When a molec¬ 
ular spin is screened by delocalized electrons, the Kondo 
effect may be observed, which gives rise to a characteris¬ 
tic anomaly in the excitation spectrum close to the Fermi 
energy [l|-l6[. Several approaches are being explored to 
control molecular Kondo systems at surfaces. These in¬ 
clude the application of electric fields 0 , mechanical de¬ 
formation [Sj, chemical modification of the ligand shell 
[jm, and electronic excitation followed by spin state 
trapping 


In purely metallic systems, the proximity of a metal 
tip in a scanning tunneling microscope (STM) has been 
shown to modify the Kondo effect through changes of the 
Co d band and exchange interaction |l5l4l7l|. It seems 
natural to extend this method to the spin states of metal 
complexes using the tip as an additional ligand. By mov¬ 
ing the tip, the ligand strength may be continuously mod¬ 
ified with the subsequent changes in the local environ¬ 
ment of the molecule, making this approach particularly 
appealing. The concomitant forces were observed to lead, 
in extreme cases, to the transfer of an atom or molecule 
from the sample to the tip [illii. 


Moreover, the electrostatic potential at the molecule 
may be modified. Although the electrical field in a 
nanoscopic junction is fundamentally interesting and 
drives current through the molecule little is known about 
it from experiments. Usually, it is assumed that the 
states of an adsorbed molecule are pinned with respect to 
the Fermi level of the substrate. A voltage drop between 
molecules and their substrate has been reported in cases 
where the molecules or parts of them were decoupled 
from a metal substrate by insulating layers [i® or by 
a particular geometric arrangement such as rod-shaped 
molecules adsorbed vertically [lH or molecular double 
deckers Q. The question of the voltage drop within a 
molecular junction has been addressed by nonequilibrium 


ab initio quantum transport calculations. These calcula¬ 
tions show that most of the voltage drops at the weakest 
bond of the junction miii. Calculations also indicate 
that a voltage drop between molecule and substrate may 
lead to negative differential resistance and two-level 
fluctuations [2^ of Cfio junctions as well as rectification 
in a molecular wire [29|. 

Here we present STM data and results of ab initio cal¬ 
culations for a metal complex on Au(lll). In the limit of 
weak interaction with the STM tip the complex exhibits 
an underscreened high-spin Kondo effect. When the tip 
is moved closer to the metal ion, the differential conduc¬ 
tance spectrum changes drastically. It turns out, how¬ 
ever, that this change is due to sub-Angstrom structural 
relaxations. Similar to a mechanically operated switch, 
the voltage drop over the junction moves from the tip- 
molecule gap to the molecule-substrate bond. 

Purified 5, 10, 15, 20-tetrakis(4-sulfonatophenyl)-21H, 
23H-porphine manganese (III) chloride (MnTPS-Cl) was 
dissolved in aqueous methanol in the presence of 1 vol% 
acetic acid. Mass-selected singly charged MnTPS cations 
originating from the dissolution of Cl ions were electro- 
sprayed onto clean Au(lll) surfaces in ultrahigh vacuum 
(UHV) [3^. Thus prepared samples were cooled to 77 K 
and subsequently transferred into a STM, which was op¬ 
erated at ^ 5 K in UHV. Tips were electrochemically 
etched from tungsten wire and further prepared in situ 
by indentation into the substrate. A sinusoidal modula¬ 
tion (2 mVrrns, 1-2 kHz) was added to the bias to record 
spectra of the differential conductance. 

Figure [T](a) schematically shows MnTPS, which hosts a 
Mn ion at the center of a tetraphenylporphyrin with four 
sulfonyl hydroxide (-SO3H) groups occupying para posi¬ 
tions of the phenyl rings. When adsorbed on Au(lll), 
the molecules appear quatrefoil shaped in STM to¬ 
pographs at positive sample voltage V [Fig. [11(b)], as 
typic ally observed for porphyrins and phthalocyanines 
[1,0, El- One of the molecular C2 axes preferably aligns 
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FIG. 1: (a) Structure of MnTPS (H atoms not shown). 

Dotted lines through pairs of opposite pyrrole units indicate 
C 2 axes, (b) STM topograph (4x4 nm^, sample voltage V = 
2.65 V, current / = 50 pA) of MnTPS agglomerated into an 
island on Au(lll). The data are displayed in a pseudo-three- 
dimensional fashion. The apparent height of the molecular 
lobes is 0.28 nm. Solid and dashed vectors indicate the 
[Oil] and [211] directions of Au(lll), respectively, (c) Spatial 
evolution of dl/dV spectra of MnTPS. Dots in the inset 
indicate the positions of the measurements (spectrum I above 
the molecular center, spectrum III at the center of a lobe). For 
clarity, spectra II and III are shifted by —0.05 and —0.1 nS, 
respectively. The solid red line depicts a Frota fit to spectrum 
I. 


along a (Oil) direction of Au(lll) within the agglomer¬ 
ates. However, deviations of ~ ±5° from this orienta¬ 
tion were observed. Slight variations in apparent heights 
among the lobes of a particular as well as of different 
molecules are likely due to incommensurability with the 
substrate lattice. This structural variability may indicate 
a fairly weak molecule-substrate interaction as expected 
for An. 


dl/dV spectra measured over various positions of the 
molecules are shown in Fig. Above the Mn center 

(spectrum I) a resonance is observed close to zero bias, 
which continuously vanishes as the tip is moved sideways 
onto a lobe (spectra II and III). Following previous re¬ 
ports from metals 0 and metal porphyrins we 

attribute the resonance as a signature of the Kondo ef¬ 
fect. The spectrum of the Mn center may be described 
well by a Frota line 33, (solid red line) leading to 
an estimated Kondo temperature Tk 100 K Hi. Tk 


FIG. 2: (a) Gonductance vs displacement Az of the tip 

from the tunneling condition (Az = 0 at K = 100 mV and 
I = 0.05 nA) towards contact at the center of a MnTPS 
molecule. At contact formation (arrow) the conductance is 
0.074 Go with Go = 2e^ jh being the conductance quantum. 
The conductance measured from different molecules scattered 
in the range of 0.07 ± 0.03Go. (b) dl jdV spectra acquired at 
Az — 0 (tunneling) and 0.33 nm (contact) at the center of a 
molecule. The bottom two curves are shifted by —0.45 and 
—0.35 nS, respectively, for clarity. Solid red lines show fits of 
Frota functions in the range of ±30 mV. 


measured from different molecules scattered in the range 
100(7) K. No systematic effect of the position of a 
molecule in the agglomerates on the line shape was ob¬ 
served. 

Next, the tip was brought closer to the center of a 
molecule until contact was reached. Data recorded for 
a molecule in a fee region of the Au(lll) reconstruction 
are shown in Fig. [2](a). Az is the displacement of the 
tip from its initial height defined by V = 100 mV and 
/ = 0.05 nA. The abrupt rise of the conductance at 
Az = 0.33 nm signals the transition to contact. 

dl/dV spectra were acquired by holding the tip at pre¬ 
determined heights over the molecular center [Fig. [2][b)]. 
Spectra recorded in the tunneling range exhibit similar 
line shapes (only spectra corresponding to Az = 0 are 
shown). At contact, a resonance near zero bias is still 
observed but surprisingly the line shape appears to be 
inverted on the voltage axis. Further approach of the tip 
from the point of contact does not change the shape of 
the spectrum. Approach of the tip beyond ^0.03 nm 
from the point of contact usually led to the destruc- 
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tion of the junction. Very similar Kondo temperatures 
100 K) are obtained from the fits of Frota line shapes 
to both the tunneling and contact spectra. Tunneling 
spectra {Az = 0) recorded after the contact experiment 
were unaltered [bottom curve in Fig.[2][b)], excluding the 
possibility that an irreversible change of the tip or the 
molecule occurred at contact. Isolated molecules, which 
were prepared by manipulation with the tip, exhibited 
tunneling characteristics virtually identical to that of the 
molecules in the agglomerates [Fig. Efb)]. However, at 
contact isolated molecules were usually unstable. 

At first glance, it may seem natural to attribute a 
changed Kondo line to a modified molecular spin state. 
On the other hand, it appears difficult to interpret the ob¬ 
served inversion of the line in such a scenario. We there¬ 
fore performed density functional theory (DFT) based ab 
initio electronic structure and transport calculations to 
investigate to what extent the Kondo physics is modified 
by the proximity of the tip. The calculation explicitly 
takes into account the strong dynamic correlations origi¬ 
nating from the Mn 3d-shell that give rise to the Kondo 
effect First, the MnTPS molecule was relaxed in 

the junction between the Au(lll) substrate and the STM 
tip, modeled by a small An pyramid built in the [111] di¬ 
rection, using VASP [m, [ 3 ^ with the exchange-correlation 
functional by Perdew, Burke and Ernzerhof [i^ and the 
dispersion correction due to Grimme [4l| (see details in 
the Supplemental Material [H^). Figures [^a) and [3][b) 
show the resulting structures in the tunneling regime and 
at contact, respectively. Without the tip, the Mn center 
is found to be at a height of 2.77 A above the surface. 
When the tip is brought closer to the Mn center, the 
molecule is lifted from the surface. At contact (tip at a 
height of 6 A above the surface), the displacement of Mn 
center is ^ 0.7 A [^ . 

In order to determine the dynamic single-particle 
broadening rc^(cc;) of the Mn 3d-levels due to the cou¬ 
pling to the organic framework of the molecule, the sub¬ 
strate, and the tip, DFT-based transport calculations 


were performed using the ANT. G code . To make the 
calculations computationally feasible we replaced the sul- 
fonatophenyl groups of the molecule by hydrogen atoms, 
which leaves the local environment of the Mn center un¬ 
changed. Figure [3l[c) shows rc^(cc;) of the orbital, 
which is the only orbital that couples directly to the 
5-type conduction electrons of the substrate (similar to 
MnPc 37, li^) and effectively binds to the s level of the 
An tip atom (details in the Supplemental Material 0). 
Hence, the 3d;22-orbital is essentially the only d-orbital 
that responds to the lifting of the molecule from the sub¬ 
strate. In the tunneling regime the broadening is fea¬ 
tureless around the Fermi energy Ep [Fig.[3][c), solid red 
line], while the coupling to the 5-level of the tip at contact 
leads to a broad feature at ^ 1.2 eV above Ep [Fig.[3][c), 
dashed blue line]. The concomitant increased broadening 
at Ep is partially compensated by reduced coupling to 
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FIG. 3: (a),(b) Schematics of relaxed structures in tunnel¬ 

ing and contact regimes. At contact, which corresponds to 
a tip height of 6 A above the substrate, the Mn center is 
lifted from the surface by 0.7 A with respect to the tunnel¬ 
ing regime. This vertical displacement is exaggerated in the 
figure for better perceptibility, (c) Broadening of 3d^2 orbital 
in the tunneling (solid red line) and contact (dashed blue 
line) regimes, (d) Calculated transmission functions T{uj) 
normalized to the maximum transmission Tmax in the tun¬ 
neling and contact (solid black lines) regimes along with fits 
using Frota (dashed cyan lines) and Fano (solid magenta lines) 
line shapes. Fits minimize the mean-square deviation over the 
interval of -100 to 100 meV. For clarity, a narrower range is 
shown. The transmission curve at contact is shifted by 1. (e) 
Calculated dl/dV curves in the tunneling (solid red line) and 
in the contact regime (dashed blue line), a is 0.02 and 0.8 in 
the tunneling and contact regimes, respectively. The contact 
spectrum is shifted by 10 /xS. 


the substrate. 

To capture correlation effects beyond the mean-field 
approximation, the Mn 3(i shell was augmented by a Hub¬ 
bard interaction term that takes into account density- 
density interactions and Hund’s rule coupling. We as¬ 
sume an intraorbital repulsion U = 4.5 eV, an interor¬ 
bital repulsion U' = 3.1 eV, and a Hund’s coupling 
Jh = 0.7 eV (HU. The one-crossing approximation 
(OCA) [ 4 J was used to solve the multiorbital Ander¬ 
son impurity model of the Mn 3d shell coupled to the 
rest of the system. The solution yields the spectral func¬ 
tion of the Mn 3d-shell Ad{uj) (see Supplemental Mate¬ 
rial 0). Only the Sd^^ orbital couples to the conduction 
electrons near Ep and thus shows a Kondo resonance. 
While the width of the resonance is comparable to that 
of the dljdV spectra recorded in the tunneling regime. 
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a slightly broader resonance appears at contact due to a 
small increase of and charge fluctuations of the 

orbital (see Supplemental Material 0). The total 
spin of the Mn 3d shell is found to be approximately 1.7, 
i.e., halfway between spin-3/2 and spin-2 systems. Hence 
MnTPS displays an underscreened high-spin Kondo ef¬ 
fect. 


To obtain the low-bias transport properties of the sys¬ 
tem we calculated the transmission function T{uj) = 
Tr[fT(co’)Gj^(ci;)f s(co’)Gm(<^)], where Gm(<^) is the cor¬ 
related Green’s function of the extended molecule that 
contains the Kondo resonance. Tt and Ts are operators 
describing the coupling to the tip and the substrate, re¬ 
spectively. Figure [3l^d) shows the calculated transmission 
curves (for both the tunneling and contact regimes) that 
exhibit asymmetric resonances similar in shape to those 
of the dl/dV spectra measured in the tunneling regime 
(Fig. Efb)). However, in contrast to the experiment no 
line shape inversion occurs at contact. Figure ^d) also 
shows hts of Fano (solid magenta lines) and Frota (dot¬ 
ted cyan) functions to the calculated transmissions for 
both tunneling and contact. Interestingly, the ht of the 
Frota line to our results, which are based on the OCA, 
is superior. Similar observations were previously made 
for spectral functions calculated using numerical renor¬ 
malization group ill, and quantum Monte Carlo 
approaches. Second order perturbation theory, on the 
other hand, results in a Fano line shape j48| . 


The difference in line shape between the calculated 
transmission functions and the measured dl/dV is un¬ 
derstood by considering the voltage drop across the junc¬ 
tion. The current through the molecule can be writ¬ 
ten in terms of the transmission function as /(K) = 
X dujT{(jj — aeV)^ where /i is the chemical poten¬ 

tial of the substrate and a is a dimensionless parameter 
determining the local chemical potential at the position 
of the molecule 23, li^. Similar to a voltage divider, 
a can be obtained from the broadenings Tt (Fg) due to 
the tip (substrate) as a = Since the transport 

through MnTPS is dominated by transmission through 
the 3dz2 orbital, for F^ and F^ we use the couplings of the 
Mn 3dz2 orbital to the tip and substrate, respectively. We 
estimate a = 0.02 with Tt ^ 0.01 eV and Tg ^ 0.54 eV 
in the tunneling regime, and a = 0.7 at contact (tip 6 A 
above the surface) with Tt ^ 0.47 eV and Tg ^ 0.2 eV. 
As expected, e. g.. from nonequilibrium quantum trans¬ 


port calculations [26|, Tt Tg in the tunneling range, 
a is small, the potential at the molecule is pinned to 
and the applied bias drops over the vacuum gap between 
the tip and the molecule. 


Figure n e) shows the calculated dl/dV spectra for 
a = 0.02 (tunneling) and a = 0.8 (contact). The lat¬ 
ter value is chosen because the estimated a = 0.7 gives a 
rather symmetric resonance (also see Supplemental Ma¬ 
terial [50|). Like in the experiments, the line shape of 


the spectrum at contact (dashed blue line) is inverted 
with respect to the tunneling regime (solid red line). 
Since no such inversion occurs in the underlying trans¬ 
mission curves [Fig.[3fd)], the observed inversion is solely 
due to the modihed local potential of the molecule. In 
other words, the vo ltag e drop has shifted to the molecule- 
substrate spacing [^. This shift is not only due to 
an increased molecule-tip coupling but also to a sub- 
Angstrom displacement of the molecule, which reduces 
Tg. Moreover, since a > 0.8 is required to produce an 
asymmetry comparable to that of the experimental data, 
we conclude that at least 80% of the voltage drops be¬ 
tween molecule and substrate in the contact regime. 

A second, less obvious difference between tunneling 
and contact is that the Kondo screening of the molec¬ 
ular spin is mainly mediated by the electron systems of 
the substrate and the tip, respectively. In the present 
experiments, the tip was covered with An, the material 
of the substrate and the Kondo temperature happens to 
be hardly changed. Occasionally, with some tips, the 
Kondo resonance disappeared in the contact range, pos¬ 
sibly because the material at the apex of those tips did 
not provide the required screening. 

The model of the voltage drop used above has been 
very useful in interpreting data from double-barrier sit¬ 
uations like tunneling to a molecule on a thin insulator. 
The data presented here for a fiat molecule between two 
metal electrodes, shows that the voltage drop can be very 
relevant (> 80% between substrate and molecule) in a 
case that is less expected. While for MnTPS {a close to 
one) an inversion of the line shape is observed, a more 
symmetric coupling {a ^ 0.5) would cause an artihcial 
broadening of the measured resonance. In summary, the 
electrical potential in molecular and atomic junctions ap¬ 
parently deserves more attention. 

Financial support of the Deutsche Forschungsge- 
meinschaft (DFG) via Sonderforschungsbereich 677 and 
Forschergruppe 1162 and 1346 is acknowledged. 


Supplemental Material 
Structure relaxations 

The structures of an MnTPS molecule have been re¬ 
laxed on a Au(III) substrate using VASP [38, 39], both 
in the absence and in the presence of a An tip (in form of 
a pyramid grown in the [III] direction) in the proximity 
of the central Mn atom. We have used the Perdew-Burke- 
Ernzerhof (PBE) functional [40] and included dispersion 
correction due to Grimme (DET-D3) [41]. The experi¬ 
mental situation has been modeled using a 3 layer slab 
of Au(III) extending over 10 x 10 lateral primitive unit 
cells. Eor the relaxations we used the F point in the 
Brillouin zone only. All crystal structures have been re¬ 
laxed until the forces acting on each atom were smaller 
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FIG. 4: Relaxed structures of MnTPS molecule on Au(lll) in the tunneling regime, (a) Top view of the molecule adsorbed on 
the fee hollow site. (b,c) Side view along the two N-Mn-N axes with a 30° inclination w.r.t. the Au surface. (d,e) Side view 
along the two N-Mn-N axes. For subsequent DFT+many-body calculations the MnP “core” of the molecule has been used. 


than 0.05 eV/A. The planar geometry of the molecule 
is strongly deformed on the surface, with the porphine 
ring forming a saddle structure. The phenyl groups ro¬ 
tate themselves out of the molecular plane by about 30° 
as shown in Fig. |U Energetically the fee hollow sites 
are found as the most favorable positions for the molec¬ 
ular center. In the absence of a tip, the distance of the 
Mn center from the Au surface is 2.77 A. When the tip is 
brought into the proximity of the Mn center, the molecule 
is lifted from the surface. For the apex atom at the 
heights of 6.0, 5.8, and 5.6 A above the surface, the Mn 
center is lifted by about 0.7, 0.5, and 0.3 A, respectively, 
from its initial position. 


DFT+OCA electronic structure and transport 
calculations 

We have performed ab initio electronic structure and 
transport calculations explicitly taking into account the 
strong dynamic correlations originating from the Mn 3d- 
shell that give rise to the Kondo effect, as described in 
previous work [35-37]. Using the ANT.G ab initio trans¬ 
port code [42] we have first performed DFT calculations 
of an embedded cluster containing the molecule, a finite 
part of the Au surface and the pyramidal Au-tip. Both 
the tip and the finite Au surface have been embedded into 
electrodes modeled by tight-binding Bethe lattices with 
realistic tight-binding parameters for bulk Au obtained 
from DFT calculations. In order to make the calcula¬ 
tions computationally feasible for the ANT.G code, we 
replaced the sulfonatophenyl groups of the molecule by 
hydrogen atoms, leaving the Mn center in a local envi¬ 
ronment similar to that of the original structure. 

From the converged DFT calculation we obtain the 


Kohn-Sham Green’s function of the embedded cluster, 

G%{aj) = [w + Ai - - St(w) - SsH]“t (1) 

where /i is the chemical potential, Hq is the Kohn-Sham 
(KS) Hamiltonian of the cluster which yields a mean-field 
description of the electronic structure of the cluster and 
and Es(<^) are the embedding self-energies which 
describe the coupling of the cluster to the bulk electrodes 
of the STM and substrate, respectively. From 
calculate the so-called hybridization function of the Mn 
3(i-shell, 

Ad{uj) = id + - [G^{uj)] ( 2 ) 

where and G^(cc;) are the projected KS Hamiltonian 
and KS Green’s function, respectively, for the Mn 3d 
states. The negative imaginary part of A^(cj) yields 
the (dynamic) single-particle broadening of the 

3(i-levels due to the coupling to the rest of the system 
comprising the organic rest of the molecule, substrate 
and the tip. Figure [5] shows the orbit ally resolved single¬ 
particle broadening for the Mn 3d-shell in the tun¬ 

neling regime. The hybridization functions resemble the 
ones calculated for manganese phthalocyanine (MnPc) 
on Pb(lll) and Ag(OOl) [37, 43], but differ in some im¬ 
portant details. The broadening of the 3(i;22-orbital is 
relatively featureless around the Fermi level indicating a 
coupling of that orbital predominantly to the 5-type con¬ 
duction electrons in the Au substrate. Near the Fermi 
level the dxz- and dyz-oihitsih do not couple directly to 
the substrate, but the appearance of Lorentzian peaks 
in the hybridization functions indicate indirect coupling 
of those orbitals to the substrate via organic ligands. In 
contrast to MnPc, these peaks are not directly at the 
Fermi level but about 1 eV above. The dx‘^-y ‘2 orbital 
has an extremely strong coupling to the organic rest of 
the molecule at — 2 eV which leads to quenching of that 











6 




FIG. 5: Orbital resolved dynamic broadening 
—ImAd(cj) for Mn Sd-shell in the tunneling regime. 


FIG. 6: Orbital resolved spectral function in the tunneling 
regime. 


orbital, while the d^^^-orbital does not couple at all to the 
substrate or the organic rest for low energies. 

In order to capture dynamic correlation effects origi¬ 
nating from the strongly interacting 3d-electrons of the 
Mn center, the Mn Sd-shell is augmented by a Hubbard 
type interaction term, 

2 ^ ^ di(j'd]^(j^ (3) 

ijklaa' 

where Uijki are the matrix elements of the effective 
Coulomb interaction for the Mn 3d-shell. The effective 
Coulomb interaction is much smaller than the bare one 
due to RPA screening by the conduction electrons in the 
substrate, tip and organic rest of the molecule. Here 
we take a simplified model interaction taking into ac¬ 
count only density-density interactions Uijij and Hand’s 
rule coupling Uijji. We assume an intra-orbital repulsion 
U = 4.5 eV, an inter-orbital repulsion U' = 3.1 eV and 
a Hand’s coupling Jh = 0.7 eV. Since the Coulomb in¬ 
teraction within the 3d-shell has already been taken into 
account on a mean-field level by the Kohn-Sham DFT 
Hamiltonian 77^, we subtract a double-counting correc¬ 
tion (DCC) in order to obtain the bare energy levels of 
the Mn 3d-shell: — Vdc- Here we choose the 

so-called fully-localized limit (FLL) DCC scheme [51]. 

The Mn 3d-levels the effective Coulomb interaction 
1-Lu and the hybridization function Ad{uj) completely de¬ 
fine a multi-orbital Anderson impurity model (AIM). We 
have used One-Crossing Approximation (OCA) to solve 
the AIM [44]. The solution of the AIM yields the spectral 
function of the Mn 3d-shell Ad{uj) which is shown in Fig. [6] 
for the tunneling regime. Only the -orbital shows a 
Kondo peak at the Fermi level. The other orbitals do 
not show Kondo peaks but have strong resonances not 
too far from the Fermi level. Overall the Mn 3d-shell is 
in a mixed valence situation with a total occupation of 
^5.5, although individual orbital-occupancies are rather 
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FIG. 7: Gorrelated transmission functions (a) in the tunneling 
regime with the tip at 7.8A above the surface, and (b) in the 
contact regime with tips at different heights above the surface 
as indicated. 


close to 1. The total spin of the Mn 3d-shell is found to be 
approximately 1.7, i.e. half-way between a Spin-3/2 and 
a Spin-2 systems. The spin in the da, 2 _^ 2 -orbital is com¬ 
pletely quenched by the strong coupling to the organic 
ligands. At contact the d; 22 -orbital is slightly lowered in 
energy due to the bonding to the tip atom, increasing the 
occupancy of that orbital to ^ 1.2. 

Using the Kramers-Kronig relation the real part of the 
correlated 3d Green’s function Gd{oj) has been calcu- 
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lated. From that we obtain the electronic self-energy 
via 'Ed = — [Gd]~^^ which captures the dynamic 

correlation effects induced by the effective Coulomb in¬ 
teraction Hu- The correlated electronic structure and 
transport properties of the whole system is then calcu¬ 
lated from the correlated cluster Green’s function, 

6c(c^) = [(Gc(c^))-' - SdH + Vdc]-\ (4) 

where the Ed and Vdc only act on the subspace of the 
Mn 3d-shell. From Gc{^) we calculate the correlated 
transmission function, 

T(:^) = IV[f T (w)G^ (w)f s {co)Gc (c^)] , (5) 

which according to Landauer yields the low-bias trans¬ 
port properties, i.e. the low-bias differential conductance 
and current [52]. Ft and Fs are the imaginary parts of 
the electrode self-energies Et and Ss, respectively. Fig¬ 
ure [3 shows the calculated transmission functions both 
in the tunneling and in the contact regime for different 
heights of the tip above the surface. For distances below 
5.8 A the transmission function is strongly broadened, 
i.e. the Kondo temperature increases strongly. Also the 
conductance increases strongly. However, junction insta¬ 
bility at higher current prohibited to make this observa¬ 
tion experimentally. 



Voltage (mV) 

FIG. 8: Calculated dl/dV spectra in the contact regime with 
the tip at GA above the substrate for different values of the 
parameter a. 

Figure [8] shows the calculated dl/dV spectra in the 
contact regime for different values of a, the dimensionless 
parameter determining the local chemical potential at the 
position of the molecule. For r* < r, (r* > r,), a ^ o 
(a —>■ 1) implying that the local potential at the molecule 
is pinned to the chemical potential of the substrate (the 
tip). 


* Electronic address: karan@physik.uni-kiel.de 


^ Present address: Key Laboratory for the Physics and 
Chemistry of Nanodevices, Department of Electronics, 
Peking University, Beijing 100871, P. R. China 
^ Electronic address: weismann@physik.uni-kiel.de 

[1] A. Zhao, Q. Li, L. Chen, H. Xiang, W. Wang, S. Pan, 

B. Wang, X. Xiao, J. Yang, J. Hou, et al.. Science 309, 
1542 (2005). 

[2] U. G. E. Perera, H. J. Kulik, V. lancu, L. G. G. V. 
Dias da Silva, S. E. Ulloa, N. Marzari, and S.-W. Hla, 
Phys. Rev. Lett. 105 , 106601 (2010). 

[3] A. Mugarza, R. Robles, C. Krull, R. Korytar, N. Lorente, 
and P. Gambardella, Phys. Rev. B 85, 155437 (2012). 

[4] E. Minamitani, N. Tsukahara, D. Matsunaka, Y. Kim, 
N. Takagi, and M. Kawai, Phys. Rev. Lett. 109 , 086602 
( 2012 ). 

[5] T. Komeda, H. Isshiki, J. Liu, Y.-E. Zhang, N. Lorente, 
K. Katoh, B. K. Breedlove, and M. Yamashita, Nat. 
Commun. 2 , 217 (2011). 

[6] A. DiLullo, S.-H. Chang, N. Baadji, K. Clark, J.-P. 
Klockner, M.-H. Prosenc, S. Sanvito, R. Wiesendanger, 

G. Hoffmann, and S.-W. Hla, Nano Lett. 12 , 3174 (2012). 

[7] V. Meded, A. Bagrets, K. Eink, R. Chandrasekar, 
M. Ruben, E. Evers, A. Bernand-Mantel, J. S. Selden- 
thuis, A. Beukman, and H. S. J. van der Zant, Phys. 
Rev. B 83, 245415 (2011). 

[8] J. J. Parks, A. R. Champagne, T. A. Costi, W. W. Shum, 
A. N. Pasupathy, E. Neuscamman, S. Elores-Torres, P. S. 
Cornaglia, A. A. Aligia, C. A. Balseiro, et al.. Science 
328 , 1370 (2010). 

[9] C. Wackerlin, D. Chylarecka, A. Kleibert, K. Muller, 

C. lacovita, E. Nolting, T. A. Jung, and N. Ballav, Nat. 
Commun. 1 , 1 (2010). 

[10] A. Strozecka, M. Soriano, J. 1. Pascual, and J. J. Palacios, 
Phys. Rev. Lett. 109 , 147202 (2012). 

[11] T. G. Gopakumar, H. Tang, J. Morillo, and R. Berndt, 
J. Am. Chem. Soc. 134 , 11844 (2012). 

[12] L. Liu, K. Yang, Y. Jiang, B. Song, W. Xiao, L. Li, 

H. Zhou, Y. Wang, S. Du, M. Ouyang, et al., Sci. Rep. 
3 , 1210 (2013). 

[13] T. G. Gopakumar, E. Matino, H. Naggert, A. Bannwarth, 
E. Tuczek, and R. Berndt, Angew. Chem. Int. Ed. 51 , 
6262 (2012). 

[14] T. Miyamachi, M. Gruber, V. Davesne, M. Bowen, 
S. Boukari, L. Joly, E. Scheurer, G. Rogez, T. K. Ya- 
mada, P. Ohresser, et al., Nat. Commun. 3 , 938 (2012). 

[15] N. Neel, J. Kroger, L. Limot, K. Palotas, W. A. Hofer, 
and R. Berndt, Phys. Rev. Lett. 98 , 016801 (2007). 

[16] N. Neel, J. Kroger, and R. Berndt, Phys. Rev. B 82 , 
233401 (2010). 

[17] J. Bork, Y. Zhang, L. Diekhoner, L. Borda, P. Simon, 
J. Kroha, P. Wahl, and K. Kern, Nat. Phys. 7 , 901 (2011). 

[18] D. M. Eigler, C. P. Lutz, and W. E. Rudge, Nature (Lon¬ 
don) 352 , 600 (1991). 

[19] G. Schull, T. Erederiksen, M. Brandbyge, and R. Berndt, 
Phys. Rev. Lett. 103 , 206803 (2009). 

[20] G. V. Nazin, S. W. Wu, and W. Ho, Proc. Natl. Acad. 
Sci. U.S.A. 102 , 8832 (2005). 

[21] B. Li, C. Zeng, J. Zhao, J. Yang, J. G. Hou, and Q. Zhu, 
J. Chem. Phys. 124 , 064709 (2006). 

[22] B. Warner, E. E. Hallak, H. Priiser, J. Sharp, M. Persson, 
A. J. Eisher, and C. E. Hirjibehedin, Nat. Nanotech. 10 , 
259 (2015). 

[23] S. Datta, W. Tian, S. Hong, R. Reifenberger, J. 1. Hen¬ 
derson, and C. P. Kubiak, Phys. Rev. Lett. 79 , 2530 








8 


(1997). 

[24] F. Matino, G. Schull, F. Kohler, S. Gabutti, M. Mayor, 
and R. Berndt, Proc. Natl. Acad. Sci. U.S.A. 108 , 961 
( 2011 ). 

[25] G. G. Liang, A. W. Ghosh, M. Paulsson, and S. Datta, 
Phys. Rev. B 69 , 115302 (2004). 

[26] K. H. Bevan, D. Kienle, H. Guo, and S. Datta, Phys. 
Rev. B 78 , 035303 (2008). 

[27] F. D. Novaes, M. Gobian, A. Garcia, P. Ordejon, 
H. Ueba, and N. Lorente, arXiv: 1101.3714. 

[28] S. Ulstrup, T. Frederiksen, and M. Brandbyge, Phys. 
Rev. B 86, 245417 (2012). 

[29] J. Taylor, M. Brandbyge, and K. Stokbro, Phys. Rev. 
Lett. 89 , 138301 (2002). 

[30] G. Hamann, R. Woltmann, 1. Hong, N. Hauptmann, 
S. Karan, and R. Berndt, Rev. Sci. Instrum. 82 , 033903 
( 2011 ). 

[31] M. Ternes, A. J. Heinrich, and W.-D. Schneider, J. Phys. 
Gondens. Matter 21, 053001 (2009). 

[32] Q. Li, S. Yamazaki, T. Eguchi, H. Kim, S.-J. Kahng, 
J. F. Jia, Q. K. Xue, and Y. Hasegawa, Phys. Rev. B 80, 
115431 (2009). 

[33] H. O. Frota, Phys. Rev. B. 45, 1096 (1992). 

[34] H. Priiser, M. Wenderoth, A. Weismann, and R. G. Ul- 
brich, Phys. Rev. Lett. 108 , 166604 (2012). 

[35] D. Jacob, K. Haule, and G. Kotliar, Phys. Rev. Lett. 
103 , 016803 (2009). 

[36] M. Karolak, D. Jacob, and A. 1. Lichtenstein, Phys. Rev. 
Lett. 107 , 146604 (2011). 

[37] D. Jacob, M. Soriano, and J. J. Palacios, Phys. Rev. B 
88, 134417 (2013). 

[38] G. Kresse and J. Hafner, J. Phys.: Gondens. Matter 6, 
8245 (1994). 

[39] G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 (1999). 

[40] J. P. Perdew, K. Burke, and M. Ernzerhof, Phys. Rev. 
Lett. 77 , 3865 (1996). 

[41] S. Grimme, J. Comput. Ghem. 27 , 1787 (2006). 

[42] D. Jacob and J. J. Palacios, J. Ghem. Phys. 134, 044118 


( 2011 ). 

[43] J. Kiigel, M. Karolak, J. Senkpiel, P.-J. Hsu, G. Sangio- 
vanni, and M. Bode, Nano Lett. 14 , 3895 (2014). 

[44] K. Haule, S. Kirchner, J. Kroha, and P. Wolfle, Phys. 
Rev. B 64, 155111 (2001). 

[45] H. O. Erota and L. N. Oliveira, Phys. Rev. B 33, 7871 
(1986). 

[46] H. Priiser, M. Wenderoth, P. E. Dargel, A. Weismann, 
R. Peters, T. Pruschke, and R. G. Ulbrich, Nat. Phys. 8, 
172 (2012). 

[47] R. N. Silver, J. E. Gubernatis, D. S. Sivia, and M. Jarrell, 
Phys. Rev. Lett. 65, 496 (1990). 

[48] B. Horvatic and V. Zlatic, Phys. Status Solid! (b) 99 , 
251 (1980). 

[49] S. W. Wu, G. V. Nazin, X. Ghen, X. H. Qiu, and W. Ho, 
Phys. Rev. Lett. 93 , 236802 (2004). 

[50] See Supplemental Material, which includes Refs. [51-52]. 

[51] M. T. Gzyzyk and G. A. Sawatzky, Phys. Rev. B 49 , 
14211 (1994). 

[52] R. Landauer, Philos. Mag. 21, 863 (1970). 

[53] Different definitions of the Kondo temperature exist in 
the literature. Eor the sake of comparability to the Kondo 
temperatures extracted from previous STM experiments 
we use the definition that relates Tk to the half width 
at half maximum of the symmetric resonance. This gives 
ksTK — 2.54F, where /cs is the Boltzmann constant and 
the width F of the Frota line is defined in REf. [53] . 

[54] At tip heights exceeding 6 A above the substrate the 
molecule was not lifted further. We thus define 6 A as 
the point of contact. 

[55] These values are slightly smaller than the ones used for 
MnPc [23,111. However, the results are hardly affected 
by changes of these parameters by up to 30%. 

[56] Reference [4^ pointed out this possibility. However, no 
corresponding experimental observation was made due to 
junction instability at high conductances. 


